Fault location method and device

ABSTRACT

The present invention relates to a method for calculating the distance to fault in a section of a power transmission network, which section is arranged with line terminals at both ends, each terminal including impedance relays (AA, BB), including the steps of, at the occurence of a fault; measuring the apparent impedances by each relay, measuring the relaying currents, determining the type of fault, checking whether the fault involves a fault resistance or not, and, if so, solving a quadratic equation for complex numbers, resolving the quadratic equation for the real and imaginary components, obtaining two quadratic equations for a fault distance in which the coefficients are real numbers, combining the two equations and obtaining a distance to fault, or, if not so, solving a quadratic equation for real numbers, obtaining two solutions for the distance to fault d 1 , d 2  and comparing the solutions as 0&lt;(d 1  or d 2 )&lt;1 pu, where pu is the length of transmission lines between the line terminals.

This application is a 371 of PCT/SE01/02771 filed on Dec. 14, 2001.

TECHNICAL FIELD

The present invention relates to fault location technique for a section of power transmission lines utilizing measurements of impedance relays installed at both line terminals.

BACKGROUND OF INVENTION

The following signals from impedance relays can, in general, be applied for fault location: phasors and/or amplitudes of relaying currents/voltages, which are composed for measurement of apparent impedance of fault loops, phasors of particular sequence components of phase currents and voltages, and measured apparent impedances of fault loops “seen” from the line terminals.

Utilization of measurements from impedance relays for fault location has been initiated by M. Sachdev and Agarwal in the papers “Accurate fault location estimates from digital impedance relay measurements,” Proceedings of Third International Conference on Developments in Power System Protection, London, 17-19 Apr. 1985, Conference Publication No. 249, pp 180-184 (paper [1]) and “A technique for estimating transmission line fault locations from digital impedance relay measurements,” IEEE Transaction on Power Delivery, Vol. 3, No. 1, January 1988, pp 121-129 (paper [2]). Their method uses the following measurements from impedance relays installed at both the line terminals:

-   -   for single phase-to-ground faults: apparent impedances of fault         loops, phasors of relaying currents and phasors of a zero         sequence current,     -   for the other fault types: apparent impedances of fault loops,         phasors of a positive and a negative sequence currents.

The method presented in papers [1-2] uses Cartesian description of the relations between the input data for different fault types. In consequence, in papers [1-2] a rather complex fault location algorithm is obtained. The algorithm contains 28 steps to be performed in a sequence dependent on a fault type.

In the method of the papers [1-2] the synchronization angle is required to be calculated. The calculations of the synchronization angle proposed in papers [1-2] are based on solving a quadratic equation for the unknown angle. As always for quadratic equations, two solutions are obtained. A specific one of them is taken further. Generally this provides the correct fault location in huge majority of the cases. However, there is no proof in papers [1-2] that this works correctly in complex configurations of transmission networks and different specifications of a fault. Further, the algorithm of the papers [1-2] is derived for a single line only, and subsequently can not handle parallel lines.

BRIEF DESCRIPTION OF INVENTION

The aim of the present invention is to remedy the above mentioned problems.

This is obtained by a method characterised by claim 1. Specific features of the present invention are characterised by the appending claims.

The method according to the present invention substantially differs from the method according to papers [1-2] with respect to the utilized input data, namely less input data is required than the above specified data for the method of papers [1-2]. The method according to the present invention utilizes the following information from impedance relays installed at both the line ends:

-   -   apparent impedances measured from the relaying voltages and         currents,     -   amplitudes of relaying currents.

It is to be understood that the above listed quantities are by nature not related to the synchronization angle. This is in contrast to the input data of the method according to papers [1-2] where the unknown synchronization angle is required for in calculations.

The advantages of the method according to the present invention can be summarised as follows:

-   -   The applied input data (apparent impedances of fault loop         signals and amplitudes of relaying currents) in the present         invention is characterized by no involvement of a         synchronization angle at all as these signals are by nature not         dependent on this angle.     -   The method according to the invention utilizes generalized         description of a fault and in consequence the developed fault         location algorithm is of compact form, i.e. the algorithm is the         same for all fault types and only one coefficient having a fixed         value depends on fault type,     -   A unique solution for a distance to fault is always provided,         regardless of the transmission network parameters and         specifications of the fault. Thus, the algorithm can be applied         for lines working in complex networks. For example, the         series-compensated transmission networks are considered as         complex grids—for which providing a completely unique solution         is especially desired. In this respect, the proposed fault         location algorithm can be easily adapted for such applications,         giving the unique solution.     -   The method is suitable for both a single line as well as for         parallel transmission lines. The algorithm is the same in         principle for both the applications. The difference is inside a         protective relay itself, where the fault loop current is         composed of somewhat different elements depending on whether it         is a single line or parallel lines arrangement. In the case of         parallel lines the fault loop current includes additional         elements, taking into account the mutual coupling between the         parallel lines for the zero sequence. Thus, the difference is         only with respect to the way of obtaining the input data of the         fault location algorithm (what is performed in a relay itself)         but not for the fault location algorithm itself. The algorithm         according to papers [1-2] is derived for a single line only         while the proposed method suits for both, a single and parallel         lines.

The proposed method is intended for use in post-fault analysis aimed at verification of operation of distance relays.

These and other aspects of the present invention and advantages with it will become apparent from the detailed description and from the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

In the following detailed description of the invention, reference will be made to the accompanying drawings, of which:

FIG. 1. shows a schematic diagram for fault location based on measurements of distance relays at the line terminals,

FIG. 2 shows a model of a two terminal transmission line for the positive sequence,

FIG. 3 shows a model of a two terminal transmission line for the negative sequence,

FIG. 4 shows a model of a two terminal transmission line for zero sequence,

FIG. 5 shows a positive sequence network for parallel lines,

FIG. 6 shows a negative sequence network of parallel lines,

FIG. 7 shows a zero sequence network of parallel lines,

FIG. 8 shows a general fault model for unsynchronized measurements of impedance relays installed at both the line terminals,

FIG. 9 shows a flow chart of the fault location algorithm according to the present invention,

FIG. 10 shows an example of a device for performing the method of the invention, and

FIGS. 11-13 show models for different types of faults.

DETAILED DESCRIPTION OF THE INVENTION

The fault location algorithm according to the present invention is capable locating faults in both single and parallel transmission lines. However, derivation of the algorithm is presented for parallel lines, FIG. 1, which is a more general case.

In case of a fault on a protected line a distance relay trips the associated circuit breaker. The decision is issued on the base of the performed impedance measurement. Apparent impedance “seen” by a relay can be determined with varieties of impedance algorithms. Here, it is considered that the impedance is determined by processing phasors of the fault loop signals: the voltage (V _(AA) _(—) _(p)—for a relay AA, V _(BA) _(—) _(p)—for a relay BA) and the current (I _(AA) _(—) _(p)—for a relay AA, I _(BA) _(—) _(p)—for a relay BA), which are composed according to the fault type. Table 1 gathers the fault loop signals, when using the traditional description and with applying the symmetrical components notation.

TABLE 1 Fault loop signals formed by a distance relay for different fault types Traditional description of the fault loop signals Description of the fault loop signals with use of symmetrical quantities ${\underset{\_}{V}}_{{AA}_{\_ p}} = {{{\underset{\_}{a}}_{1}{\underset{\_}{V}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{V}}_{AA2}} + {{\underset{\_}{a}}_{0}{\underset{\_}{V}}_{AA0}}}$ ${\underset{\_}{I}}_{{AA}_{\_ p}} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} +}$ Fault Fault loop voltage: Fault loop current: ${{{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{a}}_{0m}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{AB0}}}\quad$ type ${\underset{\_}{V}}_{{AA}_{\_ p}}$ ${\underset{\_}{I}}_{{AA}_{\_ p}}$ ${\underset{\_}{a}}_{1}$ ${\underset{\_}{a}}_{2}$ ${\underset{\_}{a}}_{0}$ a-g ${\underset{\_}{V}}_{AAa}$ ${\underset{\_}{I}}_{AAa} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{AB0}}$ 1 1 1 b-g ${\underset{\_}{V}}_{AAb}$ ${\underset{\_}{I}}_{AAL2} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{AB0}}$ a² a 1 c-g ${\underset{\_}{V}}_{AAc}$ ${\underset{\_}{I}}_{AAc} + {{\underset{\_}{k}}_{0}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{k}}_{0m}{\underset{\_}{I}}_{AB0}}$ a a² 1 a-b a-b-g a-b-c a-b-c-g ${\underset{\_}{V}}_{AAa} - {\underset{\_}{V}}_{AAb}$ ${\underset{\_}{I}}_{AAa} - {\underset{\_}{I}}_{AAb}$ 1 − a² 1 − a 0 b-c b-c-g ${\underset{\_}{V}}_{AAb} - {\underset{\_}{V}}_{AAc}$ ${\underset{\_}{I}}_{AAb} - {\underset{\_}{I}}_{AAc}$ a² − a a − a² 0 c-a c-a-g ${\underset{\_}{V}}_{AAc} - {\underset{\_}{V}}_{AAa}$ ${\underset{\_}{I}}_{AAc} - {\underset{\_}{I}}_{AAa}$ a − 1 a² − 1 0 $a = {\exp\left( {j\frac{2\pi}{3}} \right)}$ ${{\underset{\_}{k}}_{0} = \frac{{\underset{\_}{Z}}_{0L} - {\underset{\_}{Z}}_{1L}}{{\underset{\_}{Z}}_{1L}}},$ ${{\underset{\_}{k}}_{0m} = \frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1L}}},$ ${{\underset{\_}{a}}_{0m} = {0 - {{for}\quad a\quad{single}\quad{line}}}},$ ${\underset{\_}{a}}_{0m} = {{\underset{\_}{a}}_{0} - {{for}\quad{parallel}\quad{lines}}}$

The latter description of the signals is utilized for deriving the presented fault location algorithm. The circuit diagrams for a single line are shown in FIGS. 2-4, where FIG. 2 shows the positive sequence,

FIG. 3 shows the negative sequence and FIG. 4 shows the zero sequence, and for parallel lines in FIGS. 5-7.

FIG. 5 shows the positive network, FIG. 6 shows the negative network while FIG. 7 presents the zero sequence circuit for parallel lines. Note that there are no sources for both the negative and zero sequence components since the lines supply is assumed as three phase fully symmetrical. It is considered that the impedances for the negative sequence are equal to the respective impedances for the positive sequence.

In case of parallel lines it is further assumed that the impedances of both the parallel lines for particular sequences are identical (Z _(1LA)=Z _(1LB)=Z _(1L); Z _(0LA)=Z _(0LB)=Z _(0L)).

Using the sequence quantities for description of the fault loop signals one obtains for the protective relay AA, which protects the LINE A, FIG. 1: V _(AA) _(—) _(p) =a ₁ V _(AA1) +a ₂ V _(AA2) +a ₀ V _(AA0)  (1) $\begin{matrix} {{\underset{\_}{I}}_{AA\_ p} = {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{AA0}} + {{\underset{\_}{a}}_{0m}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{AB0}}}} & (2) \end{matrix}$ where:

V _(AA1), V _(AA2), V _(AA0)—positive, negative and zero sequence components of measured phase voltages,

I _(AA1), I _(AA2), I _(AA0)—positive, negative and zero sequence components of measured phase currents from the faulted line (LINE A),

I _(AB0)—zero sequence component of phase currents from the healthy line (LINE B),

a ₁, a ₂, a ₀, a _(0m)—complex coefficients dependent on a fault type (Table 1),

Z _(1L)—impedance of a whole line for the positive sequence,

Z _(0L)—impedance of a whole line for the zero sequence,

Z _(0m)—zero sequence impedance for the mutual coupling between the lines.

In all further considerations it is assumed that a time basis of measurements from the relay BA is taken as a reference. As the measurements of both relays (AA and BA) are not synchronous, the measurements of a relay AA are taken with including the unknown synchronization angle (δ).

Neglecting shunt capacitances of a line, the fault loop which is “seen” by the protective relay AA can be described with the following voltage equation: V _(AA) _(—) _(p) e ^(jδ) −dZ _(1L) I _(AA) _(—) _(p) e ^(jδ) −R _(F)( a _(F1) I _(F1) +a _(F2) I _(F2) +a _(F0) I _(F0))=0  (3) where:

d—distance to a fault,

δ—synchronization angle,

R_(F)—fault resistance,

I _(F1), I _(F2), I _(F0)—positive, negative and zero sequence components of the fault paths currents,

a _(F1), a _(F2), a _(F0)—positive, negative and zero sequence weighting coefficients dependent on the fault type (Table 2),

An analogous formula is obtained for the protective relay BA installed at the remote substation (note: there is no synchronization angle here as these measurements are taken as the reference):

  V _(BA) _(—) _(p)−(1−d) Z _(1L) I _(BA) _(—p) −R _(F)( a _(F1) +I _(F1) +a _(F2) I _(F2) +a _(F0) I _(F0))=0  (4)

Combining equations (3) and (4) results in: $\begin{matrix} {{{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}} + {\underset{\_}{I}}_{BA\_ p} - {{R_{F}\left( {{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}} \right)}\left\lbrack \quad{\frac{1}{{\underset{\_}{Z}}_{AA\_ p} - {d{\underset{\_}{Z}}_{1L}}} + \frac{1}{{\underset{\_}{Z}}_{BA\_ p} - {\left( {1 - d} \right){\underset{\_}{Z}}_{1L}}}} \right\rbrack}} = 0} & (5) \end{matrix}$

Note that the apparent impedances measured by the relays AA and BA, which are involved in (5), are not dependent on the synchronization angle, as they are determined: $\begin{matrix} {{\underset{\_}{Z}}_{AA\_ p} = {\frac{{\underset{\_}{V}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}}{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}} = \frac{{\underset{\_}{V}}_{AA\_ p}}{{\underset{\_}{I}}_{AA\_ p}}}} & \left( {5a} \right) \\ {{\underset{\_}{Z}}_{BA\_ p} = \frac{{\underset{\_}{V}}_{BA\_ p}}{{\underset{\_}{I}}_{BA\_ p}}} & \left( {5b} \right) \end{matrix}$

The crucial point of deriving the fault location algorithm relies in expressing the voltage drop across the fault resistance [the term: R_(F)(a _(F1) I _(F1)+a _(F2)+a _(F0) I _(F0))], which is involved in equations (3)-(5), in the following way: $\begin{matrix} {{R_{F}\left( {{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}} \right)} = {\frac{R_{F}}{{\underset{\_}{P}}_{FLT}}\left( {{{\underset{\_}{I}}_{AA\_ p}\quad{\mathbb{e}}^{j\quad\delta}} + {\underset{\_}{I}}_{BA\_ p}} \right)}} & (6) \end{matrix}$ where: P _(FLT)—complex coefficient dependent on the fault type (Table 3).

In consequence a general model for unsynchronized measurements of impedance relays is obtained as presented in FIG. 8.

The right hand side of equation (6) can be interpreted as the voltage drop across the equivalent fault impedance (R_(F)/P _(FLT)) when the sum of the relaying currents at both the line terminals (I _(AA) _(—) _(p)e^(jδ) and I _(BA) _(—) _(p)) flow through it.

The coefficient (P _(FLT)) introduced in equation (6) can thus be determined as the following ratio: $\begin{matrix} {{\underset{\_}{P}}_{FLT} = \frac{{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}} + {\underset{\_}{I}}_{BA\_ p}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}}} & (7) \end{matrix}$

Calculation of the coefficient (P _(FLT)) requires expressing the relaying currents at both the line terminals (I _(AA) _(—) _(p)e^(jδ) and I _(BA) _(—) _(p)) with use of the symmetrical quantities [equations (1), (2)—the coefficients see Table 1] and assuming values for the weighting coefficients (a _(F1), a _(F2), a _(F0)). Table 2 below gathers three sets of these weighting coefficients excluding the zero sequence component of the fault path current (a _(F0)=0). They differ in preference of using a particular component (the 1-st set prefers the negative sequence, the 2-nd prefers the positive sequence, the 3-rd has no preference).

Moreover, for determining the coefficient (P _(FLT)) it is utilized (with neglecting the line shunt capacitances) that:

-   -   for any of the sequence (positive, negative and the zero         sequence, respectively) the fault path current is received as a         sum of the currents from both the line terminals (the first         three equations of (8)) and     -   the zero sequence components of the healthy line currents at         both the stations (A and B) have opposite directions (the fourth         equation of (8)):         I _(F1) =I _(AA1) e ^(jδ) +I _(BA1)         I _(F2) =I _(AA2) e ^(jδ) +I _(BA2)         I _(F0) =I _(AA0) e ^(jδ) +I _(BA0)         I _(AB0) e ^(jδ) +I _(BB0)=0  (8)

Finally, utilizing formulae (1)-(2) together with Table 1, the coefficient (P _(FLT)) can be calculated with the relations of equation (8) and Table 2. Regardless of which of the sets (1-st, 2-nd or 3-rd) of weighting coefficients from Table 2 is used the coefficient (P _(FLT)) is of a fixed value (complex or real number) dependent on a fault type, Table 3 below. Examples of determining the coefficient (P _(FLT)) are shown in the Appendix 1.

TABLE 2 Sets of the weighting coefficients (a_(F1), a_(F2), a_(F0)) for different fault types Voltage drop across the fault path: V_(F) = R_(F)(a_(F1)I_(F1) + a_(F2)I_(F2) + a_(F0)I_(F0)) Fault 1-st set 2-nd set 3-rd set type a_(F1) a_(F2) a_(F0) a_(F1) a_(F2) a_(F0) a_(F1) a_(F2) a_(F0) a-g 0 3 0 3 0 0 1.5 1.5 0 b-g 0 3a 0 3a² 0 0 1.5a² 1.5a 0 c-g 0 3a² 0 3a 0 0 1.5a 1.5a² 0 a-b 0 1 − a 0 1 − a² 0 0 0.5(1 − 0.5(1 − 0 a²) a) b-c 0 a − 0 a² − a 0 0 0.5(a² − 0.5(a − 0 a² a) a²) c-a 0 a² − 1 0 a − 1 0 0 0.5(a − 0.5(a² − 0 1) 1) a-b-g 1 − 1 − a 0 1 − a² 1 − a 0 1 − a² 1 − a 0 a² b-c-g a² − a − 0 a² − a a − 0 a² − a a − a² 0 a a² a² c-a-g a − a² − 1 0 a − 1 a² − 1 0 a − 1 a² − 1 0 1 a-b-c-g 1 − 0 0 1 − a² 0 0 1 − a² 0 0 a² (a-b-c)

TABLE 3 Values of the coefficient P _(FLT) for different fault types a-b-g b-c-g a-g a-b c-a-g FAULT b-g b-c a-b-c TYPE c-g c-a a-b-c-g ${\underset{\_}{P}}_{FLT}$ $\frac{{2{\underset{\_}{Z}}_{1L}} + {\underset{\_}{Z}}_{0L}}{3{\underset{\_}{Z}}_{1L}}$ 2 1

Introducing the coefficient (P _(FLT)) into (5) one obtains the following quadratic equation for complex numbers: $\begin{matrix} {{{{\underset{\_}{D}}_{2}d^{2}} + {{\underset{\_}{D}}_{1}d} + {\underset{\_}{D}}_{0} + {\frac{R_{F}}{{\underset{\_}{P}}_{FLT}}{\underset{\_}{N}}_{0}}} = 0} & (9) \end{matrix}$ where: D ₂=Z _(1L) Z _(1L) D ₁ =Z _(1L) Z _(BA) _(—) _(p) −Z _(1L) Z _(1L) −Z _(1L) Z _(AA) _(—) _(p) D ₀ =Z _(1L) Z _(AA) _(—) _(p) −Z _(AA) _(—) _(p) Z _(BA) _(—) _(p) N ₀ =Z _(AA) _(—) _(p) +Z _(BA) _(—) _(p) −Z _(1L)

There are two unknowns in (9): d—distance to a fault, R_(F)—fault resistance. Solution of (9) depends on whether a fault is solid or resistive.

For ideally solid faults it is satisfied: N ₀ =Z _(AA) _(—) _(p) +Z _(BA) _(—) _(p) −Z _(1L)=0  (10)

It is to be observed that the condition (10) is satisfied also for any (solid or resistive) fault occurring outside a line. However, for such the outside faults the directional elements of distance relays are blocking the operation of relays.

Identification of solid faults or faults involving resistance of small value (a fraction of Ω) can be performed on the base of a condition that a right-hand side of (10) assumes also small values (in comparison to amplitude of a positive sequence impedance of a line): | Z _(AA) _(—) _(p) +Z _(BA) _(—) _(p) −Z _(1L) |<γ|Z _(1L)|  (11) where:

γ—coefficient being a small fraction (for example assuming γ=0.005 will provide accuracy of fault location of the range of 0.5%).

Thus, in the case of satisfying equation (11) a distance to a fault can be determined by solving the following quadratic equation, obtained from equation (9) by neglecting the term involving a fault resistance: D ₂ d ² +D ₁ d+D ₀=0  (12)

Resolving equation (12) for the real and imaginary components one obtains two quadratic equations for a fault distance in which the coefficients are real numbers. Combining these two equations the solution of equation (12) for a distance to fault, d_(s), is obtained as: $\begin{matrix} {d_{s} = \frac{{{- {{real}\left( {\underset{\_}{D}}_{0} \right)}}{{imag}\left( {\underset{\_}{D}}_{2} \right)}} + {{{imag}\left( {\underset{\_}{D}}_{0} \right)}{{real}\left( {\underset{\_}{D}}_{2} \right)}}}{{{{real}\left( {\underset{\_}{D}}_{1} \right)}{{imag}\left( {\underset{\_}{D}}_{2} \right)}} - {{{imag}\left( {\underset{\_}{D}}_{1} \right)}{{real}\left( {\underset{\_}{D}}_{2} \right)}}}} & \left( {12a} \right) \end{matrix}$ where: the complex coefficients (D ₂, D ₁, D ₀) are defined in equation (9).

In contrary, if equation (11) is not satisfied, a fault resistance has to be taken into account for fault distance estimation. The equation (9) is for complex numbers and contains two unknowns:

d—a distance to a fault, R_(F)—fault resistance. Resolving equation (9) into two equations for real and imaginary components, respectively, results in the following quadratic equation for a sought fault distance: F _(A)(d)=A ₂ d ² +A ₁ d+A ₀=0  (13) where: $\begin{matrix} {{A_{2} = {{{{real}\left( {\underset{\_}{D}}_{2} \right)}{{imag}\left( \frac{{\underset{\_}{N}}_{0}}{{\underset{\_}{P}}_{FLT}} \right)}} - {{{imag}\left( {\underset{\_}{D}}_{2} \right)}{{real}\left( \frac{{\underset{\_}{N}}_{0}}{{\underset{\_}{P}}_{FLT}} \right)}}}}{A_{1} = {{{{real}\left( {\underset{\_}{D}}_{1} \right)}{{imag}\left( \frac{{\underset{\_}{N}}_{0}}{{\underset{\_}{P}}_{FLT}} \right)}} - {{{imag}\left( {\underset{\_}{D}}_{1} \right)}{{real}\left( \frac{{\underset{\_}{N}}_{0}}{{\underset{\_}{P}}_{FLT}} \right)}}}}{A_{0} = {{{{real}\left( {\underset{\_}{D}}_{0} \right)}{{imag}\left( \frac{{\underset{\_}{N}}_{0}}{{\underset{\_}{P}}_{FLT}} \right)}} - {{{imag}\left( {\underset{\_}{D}}_{0} \right)}{{real}\left( \frac{{\underset{\_}{N}}_{0}}{{\underset{\_}{P}}_{FLT}} \right)}}}}} & \left( {13a} \right) \end{matrix}$

Equation (13), with taking equation (13a), gives two solutions for a distance to fault (d₁, d₂): $\begin{matrix} {{d_{1} = \frac{{- A_{1}} - \sqrt{A_{1}^{2} - {4A_{2}A_{0}}}}{2\quad A_{2}}}{d_{2} = \frac{{- A_{1}} + \sqrt{A_{1}^{2} - {4A_{2}A_{0}}}}{2A_{2}}}} & \left( {13b} \right) \end{matrix}$

If only one single solution satisfies equation (13b), then that solution indicates a fault in a line: 0<(d₁ or d₂)<1 pu  (13c) then the solution satisfying (13c) is taken as the valid, while the other solution (indicating a fault outside a line) is rejected.

Thus, a fault location is performed in such cases (13c) by using only apparent impedances measured by distance relays from both the line terminals. It should be appreciated that these impedances are by nature not dependent on the synchronization angle.

In contrast, if both the solutions of equation (13b) are within a line length:

 0<(d₁ and d₂)<1 pu  (13d)

then a problem of selecting the valid result arises.

According to the present invention such a selection is performed by utilizing information contained in amplitudes of relaying currents from both the line terminals. The ratio of amplitudes of relaying currents (Q) is related with a distance to fault (d) and measured apparent impedances (Z _(AA) _(—) _(p), Z _(BA) _(—) _(p)) as follows: $\begin{matrix} {Q = {{\frac{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}}{{\underset{\_}{I}}_{BA\_ p}}} = {\frac{{\underset{\_}{I}}_{AA\_ p}}{{\underset{\_}{I}}_{BA\_ p}} = {\frac{{\underset{\_}{Z}}_{BA\_ p} - {\left( {1 - d} \right){\underset{\_}{Z}}_{1L}}}{{\underset{\_}{Z}}_{AA\_ p} - {d{\underset{\_}{Z}}_{IL}}}}}}} & (14) \end{matrix}$

After some rearrangements the equation (14) transforms to the following quadratic equation (derivation is provided in APPENDIX 2): F _(B)(d)=B ₂ d ² +B ₁ d+B ₀=0  (15) where: B ₂ =|Z _(1L)|²(1−Q ²) B ₁=2real{( Z _(BA) _(—) _(p) −Z _(1L)) Z _(1L)*}+2Q ²real{ Z _(AA) _(—) _(p) Z _(1L)*} B ₀ =|Z _(BA) _(—) _(p) −Z _(1L)|² −Q ²|Z _(AA) _(—) _(p)|²  (15a)

| | and *—denote absolute value and conjugate of a complex number, respectively.

Equation (15), taking into account equation (15a), gives two solutions for a distance to fault (d₃, d₄): $\begin{matrix} {{d_{3} = \frac{{- B_{1}} - \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}}{d_{4} = \frac{{- B_{1}} + \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}}} & \left( {15b} \right) \end{matrix}$

The solutions, taken out of all four solutions: (d₁, d₂)—eq. (13b), (d₃, d₄)—eq. (15b), which coincide (d_(i)−d_(j)=0 where: i=1 or 2, j=3 or 4) give the valid solution for the fault distance (d_(v)). In practice, there are some errors in fault distance estimation (especially due to shunt capacitances of a line, which are not included at this stage of the algorithm). Thus, instead of the ideal condition (d_(i)−d_(j)=0) the following can be applied: if |d _(i) −d _(j)|=min then $\begin{matrix} {d_{v} = \frac{d_{i} + d_{j}}{2}} & (16) \end{matrix}$ where: i=1 or 2, j=3 or 4.

Note: the valid result (d_(v)) is taken as an average of the solutions (d_(i), d_(j)), which coincide.

FIG. 10 shows an embodiment of a device for determining the distance from a station, at one end of a transmission line, until the occurrence of a fault on the transmission line according to the described method. It comprises relays 12,14 arranged at both ends of the transmission line section, where the relays shown protect line A. An equivalent set of relays is used but not shown in FIG. 10 for simplicity. Each relay comprises measuring devices 1 to 3 for continuous measurement of the voltage and currents from both lines A and B. Each relay is connected to breakers 18 for each line and is equipped with means for detecting a fault and breaking the lines. Depending on the type of relay used, it may or may not comprise further devices. If the relay is of a conventional type, it may only be able to break the lines. The measured values of the currents and voltages are in that case passed to a calculating unit 30 comprising the fault location algorithm according to the invention. The measured signals are passed through measurement converters 4 to 6. The impedance may as mentioned earlier either be determined by the calculating unit by processing the phasors. The impedance may also be determined by some other means outside the calculating unit.

For a more modern type of relay, these may also comprise measurement converters, as indicated in FIG. 10, and means for classification of the fault type. Each relay comprises electronic equipment for processing phasors of the fault loop signals whereby the impedance is determined. Each relay 12, 14 is connected to the fault location calculating unit 30 via electrical connections, and the relays provide information regarding the fault type ft, the impedance Z and the relaying currents I.

The calculating unit is provided with the calculation algorithms described, programmed for the processes needed for calculating the distance to fault and the fault resistance. The calculating unit may as mentioned above also be provided with known values such as the impedance of the line. When the calculating unit has determined the distance to fault, it is displayed on the device and/or sent to remotely located display means. A printout of the result may also be provided. In addition to signalling the fault distance, the device can produce reports, in which are recorded measured values of the currents of both lines, voltages, type of fault and other associated with a given fault at a distance.

The information in the form of a result for d_(v) from the fault location system may also be embodied as a data signal for communication via a network to provide a basis for control action. The distance d_(v) may be sent as a signal as a control action such as: automatic notification to operational network centres of fault and it's location or to automatically start calculations to determine journey time to location, which repair crew shall be dispatched to site, possible time taken to execute a repair, calculate which vehicles or crew member may be needed, how many shifts work per crew will be required and the like actions.

The calculating unit may comprise filters for filtering the signals, A/D-converters for converting and sampling the signals and a micro processor. The micro processor comprises a central processing unit CPU performing the following functions: collection of measured values, processing of measured values, calculation of distance to fault and output of result from calculation. The micro processor further comprises a data memory and a program memory.

A computer program for carrying out the method according to the present invention is stored in the program memory. It is to be understood that the computer program may also be run on general purpose computer instead of a specially adapted computer.

The software includes computer program code elements or software code portions that make the computer perform the said method using the equations, algorithms, data and calculations previously described. It may also be run in a distributed way over a network. A part of the program may be stored in a processor as above, but also in a RAM, ROM, PROM or EPROM chip or similar. The program in part or in whole may also be stored on, or in, other suitable computer readable medium such as a magnetic disk, CD-ROM or DVD disk, hard disk, magneto-optical memory storage means, in volatile memory, in flash memory, as firmware, or stored on a data server.

It is to be noted that the embodiment of the invention described and shown in the drawings is to be regarded as a non-limiting example of the invention and that the scope of protection is defined by the patent claims.

Appendix 1—Derivation of the Coefficient (P _(FLT))

The coefficient (P _(FLT)) is defined in (7) as: $\begin{matrix} {{\underset{\_}{P}}_{FLT} = \frac{{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}} + {\underset{\_}{I}}_{BA\_ p}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}}} & (7) \end{matrix}$

The derivation for a single phase-to-ground (a-g) fault, FIG. 10, and for a phase-to-phase fault (a-b), FIG. 11, is performed as follows.

-   -   a-g fault:

Derivation of the weighting coefficients for particular sequence components

(the 1-st set in Table 2):

Taking into account that in the healthy phases: I _(Fb)=I _(Fc)=0 gives: ${\underset{\_}{I}}_{F1} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a{\underset{\_}{I}}_{Fb}} + {a^{2}{\underset{\_}{I}}_{Fc}}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a0} + {a^{2}0}} \right)} = {\frac{1}{3}{\underset{\_}{I}}_{Fa}}}}$ ${\underset{\_}{I}}_{F2} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a^{2}{\underset{\_}{I}}_{Fb}} + {a{\underset{\_}{I}}_{Fc}}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a^{2}0} + {a0}} \right)} = {\frac{1}{3}{\underset{\_}{I}}_{Fa}}}}$ ${\underset{\_}{I}}_{F0} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {\underset{\_}{I}}_{Fb} + {\underset{\_}{I}}_{Fc}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + 0 + 0} \right)} = {\frac{1}{3}{\underset{\_}{I}}_{Fa}}}}$

The sequence components are related: I _(F1)=I _(F2)=I _(F0)

and finally: I _(F)=I _(Fa)=3I _(F2),

thus a _(F1)=0, a _(F2)=3, a _(F0)=0 (as in Table 2)

Derivation of the coefficient P _(FLT) for a-g fault: $\begin{matrix} {{\underset{\_}{P}}_{FLT} = \frac{{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}} + {\underset{\_}{I}}_{BA\_ p}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}}} \\ {= {{by}\quad{using}\quad{the}\quad{formulas}\quad(1)\quad{and}\quad(2)\quad{we}\quad{continue}\quad{as}\quad{follows}}} \\ {= {\frac{\begin{matrix} {{\left( {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{AA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{AA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{AA0}}} \right){\mathbb{e}}^{j\quad\delta}} +} \\ \left( {{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{BA1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{BA2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{BA0}}} \right) \end{matrix}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}} +}} \\ {\frac{{{\underset{\_}{a}}_{0m}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{AB0}{\mathbb{e}}^{j\quad\delta}} + {{\underset{\_}{a}}_{0m}\frac{{\underset{\_}{Z}}_{0m}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{BB0}}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}}} \\ {= {{by}\quad{using}\quad{the}\quad{formula}\quad(8)\quad{we}\quad{continue}\quad{as}\quad{follows}}} \\ {= \frac{{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{0}\frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{F0}}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}}} \\ {= {{by}\quad{taking}\quad{the}\quad{coefficients}\quad{from}\quad{Tables}\quad 1\quad{and}\quad 2\quad{we}}} \\ {{continue}\quad{as}\quad{follows}} \\ {= \frac{{\underset{\_}{I}}_{F1} + {\underset{\_}{I}}_{F2} + {\frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}{\underset{\_}{I}}_{F0}}}{{0{\underset{\_}{I}}_{F1}} + {3{\underset{\_}{I}}_{F2}} + {0{\underset{\_}{I}}_{F0}}}} \\ {= \frac{\left( {1 + 1 + \frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}} \right){\underset{\_}{I}}_{F1}}{3{\underset{\_}{I}}_{F2}}} \\ {= \frac{\left( {1 + 1 + \frac{{\underset{\_}{Z}}_{0L}}{{\underset{\_}{Z}}_{1L}}} \right)}{3}} \\ {= \frac{{2\quad{\underset{\_}{Z}}_{1L}} + {\underset{\_}{Z}}_{0L}}{3{\underset{\_}{Z}}_{1L}}} \end{matrix}$

-   -   as in Table 3.     -   a-b fault:

Derivation of the weighting coefficients for particular sequence components

(the 1-st set in Table 2):

The fault current can be expressed as: I _(F)=I _(Fa) or as: ${\underset{\_}{I}}_{F} = {\frac{1}{2}\left( {{\underset{\_}{I}}_{Fa} - {\underset{\_}{I}}_{Fb}} \right)}$

Taking into account that in the healthy phase: I _(Fc)=0 and for the faulted phases: I _(Fb)=−I _(Fa), giving: $\begin{matrix} {{\underset{\_}{I}}_{F1} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a{\underset{\_}{I}}_{Fb}} + {a^{2}{\underset{\_}{I}}_{Fc}}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a\left( {- {\underset{\_}{I}}_{Fa}} \right)} + {a^{2}0}} \right)} = {\frac{1}{3}\left( {1 - a} \right){\underset{\_}{I}}_{Fa}}}}} \\ {{\underset{\_}{I}}_{F2} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a^{2}{\underset{\_}{I}}_{Fb}} + {a{\underset{\_}{I}}_{Fc}}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {a^{2}\left( {- {\underset{\_}{I}}_{Fa}} \right)} + {a0}} \right)} = {\frac{1}{3}\left( {1 - a^{2}} \right){\underset{\_}{I}}_{Fa}}}}} \\ {{\underset{\_}{I}}_{F0} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + {\underset{\_}{I}}_{Fb} + I_{Fc}} \right)} = {{\frac{1}{3}\left( {{\underset{\_}{I}}_{Fa} + \left( {- {\underset{\_}{I}}_{Fa}} \right) + 0} \right)} = 0}}} \end{matrix}$

The relation between I _(F1) and I _(F2) is thus: $\frac{{\underset{\_}{I}}_{F1}}{{\underset{\_}{I}}_{F2}} = {\frac{\frac{1}{3}\left( {1 - a} \right){\underset{\_}{I}}_{Fa}}{\frac{1}{3}\left( {1 - a^{2}} \right){\underset{\_}{I}}_{Fa}} = \frac{\left( {1 - a} \right)}{\left( {1 - a^{2}} \right)}}$

Finally: ${{\underset{\_}{I}}_{F} = {{\underset{\_}{I}}_{Fa} = {{\frac{3}{\left( {1 - a^{2}} \right)}{\underset{\_}{I}}_{F2}} = {\left( {1 - a} \right){\underset{\_}{I}}_{F2}}}}},$ thus: a _(F1)=0, a _(F2)=1−a, a _(F0)=0 (as in Table 2)

Derivation of the coefficient P _(FLT) for a-b fault $\begin{matrix} {{\underset{\_}{P}}_{FLT} = \frac{{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\quad\delta}} + {\underset{\_}{I}}_{BA\_ p}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}} + {{\underset{\_}{a}}_{F0}{\underset{\_}{I}}_{F0}}}} \\ {= {{as}\quad{for}\quad{this}\quad{fault}\quad{there}\quad{is}\quad{no}\quad{zero}\quad{sequence}\quad{component}}} \\ {= \frac{{{\underset{\_}{a}}_{1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{2}{\underset{\_}{I}}_{F2}}}{{{\underset{\_}{a}}_{F1}{\underset{\_}{I}}_{F1}} + {{\underset{\_}{a}}_{F2}{\underset{\_}{I}}_{F2}}}} \\ {= {{by}\quad{taking}\quad{the}\quad{coefficients}\quad{from}\quad{Tables}\quad 1\quad{and}\quad 2\quad{we}}} \\ {continue} \\ {= \frac{{\left( {1 - a^{2}} \right){\underset{\_}{I}}_{F1}} + {\left( {1 - a} \right){\underset{\_}{I}}_{F2}}}{{0{\underset{\_}{I}}_{F1}} + {\left( {1 - a} \right){\underset{\_}{I}}_{F2}}}} \\ {= \frac{{\left( {1 - a^{2}} \right)\frac{{\underset{\_}{I}}_{F1}}{{\underset{\_}{I}}_{F2}}} + \left( {1 - a} \right)}{\left( {1 - a} \right)}} \\ {= {{by}\quad{utilizing}\quad{the}\quad{relation}\quad{between}\quad{\underset{\_}{I}}_{F1}\quad{and}\quad{\underset{\_}{I}}_{F2}\quad{we}\quad{obtain}}} \\ {= {\frac{{\left( {1 - a^{2}} \right)\frac{\left( {1 - a} \right)}{\left( {1 - a^{2}} \right)}} + \left( {1 - a} \right)}{\left( {1 - a} \right)} = 2}} \end{matrix}$

-   -   as in Table 3.

Values of the coefficient (P _(FLT)) for phase-to-phase-to-ground faults and for three phase faults (Table 3) were calculated by assuming the models of these faults as in examples shown in FIG. 12.

Appendix 2—Derivation of the Coefficients in Equation (15a)

The coefficients from equation (15a), which are introduced in the quadratic equation (15) are derived as follows.

Derivation starts from eq. (14): $\begin{matrix} {Q = {{\frac{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\delta}}{{\underset{\_}{I}}_{BA\_ p}}} = {{\frac{{\underset{\_}{I}}_{AA\_ p}}{{\underset{\_}{I}}_{BA\_ p}}} = {\frac{{\underset{\_}{Z}}_{BA\_ p} - {\left( {1 - d} \right){\underset{\_}{Z}}_{IL}}}{{\underset{\_}{Z}}_{AA\_ p} - {d{\underset{\_}{Z}}_{IL}}}}}}} & (14) \end{matrix}$

Introducing: Z _(AA) _(—) _(p) =R _(AA) _(—) _(p) +jX _(AA) _(—) _(p) Z _(BA) _(—) _(p) =R _(BA) _(—) _(p) +jX _(BA) _(—) _(p) Z _(1L) =R _(1L) +jX ₁ _(—) _(L)  (14a) eq. (14) can be written down as: $\begin{matrix} {Q = {\frac{\left( {R_{BA\_ p} - R_{IL} + {dR}_{IL}} \right) + {j\left( {X_{BA\_ p} - X_{IL} + {dX}_{IL}} \right)}}{\left( {R_{AA\_ p} - {dR}_{IL}} \right) + {j\left( {X_{AA\_ p} - {dX}_{IL}} \right)}}}} & \left( {14b} \right) \end{matrix}$

Calculating absolute value of the right-hand side of equation (14b) and rising both the sides to the second power one obtains: $\begin{matrix} {Q^{2} = \frac{{\left\{ {R_{IL}^{2} + X_{IL}^{2}} \right\} d^{2}} + {\left\{ {{2{R_{IL}\left( {R_{BA\_ p} - R_{IL}} \right)}} + {2{X_{IL}\left( {X_{BA\_ p} - X_{IL}} \right)}}} \right\} d} + \left\{ {\left( {R_{BA\_ p} - R_{IL}} \right)^{2} + \left( {X_{BA\_ p} - X_{IL}} \right)^{2}} \right\}}{{\left\{ {R_{IL}^{2} + X_{IL}^{2}} \right\} d^{2}} - {\left\{ {{2R_{IL}R_{AA\_ p}} + {2X_{IL}X_{AA\_ p}}} \right\} d} + \left\{ {R_{AA\_ p}^{2} + X_{AA\_ p}^{2}} \right\}}} & \left( {14c} \right) \end{matrix}$

Utilizing the following relations, which are valid for any complex numbers: A=A ₁ +jA ₂ , B=B ₁ +jB ₂: {A ₁}² +{A ₂}² =|A| ²  (14d) {B ₁}² +{B ₂}² =|B| ²  (14e) real{ AB *}=real{(A ₁ +jA ₂)(B ₁ −jB ₂)}=real{(A ₁ B ₁ +A ₂ B ₂)+(A ₂ B ₁ −A ₁ B ₂)}=A ₁ B ₁ +A ₂ B ₂  (14f) where:

X* denotes conjugate of X,

|X| denotes absolute value of X,

eq. (14c) yields to the equation (15): F _(B)(d)=B ₂ d ² +B ₁ d+B ₀=0  (15) where: B ₂ =|Z _(1L)|²(1−Q ²) B ₁=2real{( Z _(BA) _(—) _(p) −Z _(1L)) Z _(1L)*}+2Q ²real{ Z _(AA) _(—) _(p) Z _(1L)*}  (15a) B ₀ =|Z _(BA) _(—) _(p) −Z _(1L)|² −Q ² |Z _(AA) _(—) _(p)|² 

1. Method for calculating the distance to fault in a section of a power transmission network (AA-BB), which section is arranged with line terminals at both ends, each terminal comprising impedance relays (AA, BB), comprising the steps of, at the occurrence of a fault: measuring the apparent impedances by each relay (Z _(AA) _(—) _(p), Z _(BA) _(—) _(p)), measuring the relaying currents (|I _(AA) _(—) _(p)|, |I _(BA) _(—) _(p)|), determining the type of fault, characterised in checking whether the fault involves a fault resistance or not according to |Z _(AA) _(—) _(p)+Z _(BA) _(—) _(p)−Z _(1L)|<γ|Z _(1L)| where Z _(1L)—impedance of a whole line for the positive sequence γ—small fraction coefficient <0,01, and, if fulfilled, solving a quadratic equation for complex numbers according to D ₂ d ² +D ₁ d+D ₀=0 where D ₂, D ₁, D ₀ are complex coefficients, resolving the quadratic equation for the real and imaginary components obtaining two quadratic equations for a fault distance in which the coefficients are real numbers, combining the two equations and obtaining a distance to fault as: $d_{s} = \frac{{{- {{real}\left( {\underset{\_}{D}}_{0} \right)}}{{imag}\left( {\underset{\_}{D}}_{2} \right)}} + {{{imag}\left( {\underset{\_}{D}}_{0} \right)}{{real}\left( {\underset{\_}{D}}_{2} \right)}}}{{{real}\left( {\underset{\_}{D}}_{1} \right){{imag}\left( {\underset{\_}{D}}_{2} \right)}} - {{{imag}\left( {\underset{\_}{D}}_{1} \right)}{{real}\left( {\underset{\_}{D}}_{2} \right)}}}$ or, if not fulfilled solving a quadratic equation for real numbers according to F _(A)(d)=A ₂ d ² +A ₁ d+A ₀=0 obtaining two solutions for the distance to fault according to $\begin{matrix} {d_{1} = \frac{{- A_{1}} - \sqrt{A_{1}^{2} - {4A_{2}A_{0}}}}{2A_{2}}} \\ {d_{2} = \frac{{- A_{1}} + \sqrt{A_{1}^{2} - {4A_{2}A_{0}}}}{2A_{2}}} \end{matrix}$ comparing the solutions as 0<(d ₁ or d ₂)<1pu where pu is the length of transmission lines between the line terminals.
 2. Method according to claim 1, characterised in the further steps, if the comparison 0<(d₁ and d₂)<1 pu is fulfilled: utilizing information contained in the amplitudes of relaying currents from both line terminals according to ${Q = {{\frac{{\underset{\_}{I}}_{AA\_ p}{\mathbb{e}}^{j\delta}}{{\underset{\_}{I}}_{BA\_ p}}} = {{\frac{{\underset{\_}{I}}_{AA\_ p}}{{\underset{\_}{I}}_{BA\_ p}}} = {\frac{{\underset{\_}{Z}}_{BA\_ p} - {\left( {1 - d} \right){\underset{\_}{Z}}_{IL}}}{{\underset{\_}{Z}}_{AA\_ p} - {d{\underset{\_}{Z}}_{IL}}}}}}},$ solving a quadratic equation according to F _(B)(d)=B ₂ d ² +B ₁ d+B ₀=0 where: B ₂ =|Z _(1L)|²(1−Q ²) B ₁=2real{(Z _(BA) _(—) _(p) −Z _(1L)) Z _(1L)*}+2Q ²real{Z _(AA) _(—) _(p) Z _(1L)*} B ₀ =|Z _(BA) _(—) _(p) −Z _(1L)|² −Q ² |Z _(AA) _(—) _(p)|², where | | and *—denote absolute value and conjugate of a complex number, respectively, obtaining two two solutions for a distance to fault (d₃, d₄): $\begin{matrix} {d_{3} = \frac{{- B_{1}} - \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}} \\ {d_{4} = \frac{{- B_{1}} + \sqrt{B_{1}^{2} - {4B_{2}B_{0}}}}{2B_{2}}} \end{matrix}$ comparing all four solutions: (d₁, d₂, d₃, d₄) which coincide with (d_(i)−d_(j)=0 where: i=1 or 2, j=3 or 4), thereby obtaining the valid solution for the fault distance (d_(v)).
 3. Method according to claim 2, characterised in the further step for taking into consideration of shunt capacitances of a line, comprising performing: if |d _(i) −d _(j)|=min then $d_{v} = \frac{d_{i} + d_{j}}{2}$ where: i=1 or 2, j=3 or
 4. 4. Device for calculating the distance to fault in a section of a power transmission network, which section is arranged with line terminals at both ends, each terminal comprising impedance relays (AA, BB), comprising: means for measuring the apparent impedances by each relay (Z _(AA) _(—) _(p), Z _(BA) _(—) _(p)), means for measuring the relaying currents (|I _(AA) _(—) _(p)|, |I _(BA) _(—) _(p)), means determining the type of fault, means for checking whether the fault involves a fault resistance or not according to |Z _(AA) _(—) _(p)+Z _(BA) _(—) _(p)−Z _(1L)|<γ|Z _(1L)| where Z _(1L)—impedance of a whole line for the positive sequence, /γ—small fraction coefficient<0,01, and, if fulfilled, means for solving a quadratic equation for complex numbers according to D ₂ d ² +D ₁ d+D ₀=0 where D ₂, D ₁, D ₀ are complex coefficients, means for resolving the quadratic equation for the real and imaginary components means for obtaining two quadratic equations for a fault distance in which the coefficients are real numbers, means for combining the two equations and obtaining a distance to fault as: $d_{s} = \frac{{{- {{real}\left( {\underset{\_}{D}}_{0} \right)}}{{imag}\left( {\underset{\_}{D}}_{2} \right)}} + {{{imag}\left( {\underset{\_}{D}}_{0} \right)}{{real}\left( {\underset{\_}{D}}_{2} \right)}}}{{{real}\left( {\underset{\_}{D}}_{1} \right){{imag}\left( {\underset{\_}{D}}_{2} \right)}} - {{{imag}\left( {\underset{\_}{D}}_{1} \right)}{{real}\left( {\underset{\_}{D}}_{2} \right)}}}$ or, if not fulfilled means for solving a quadratic equation for real numbers according to F _(A)(d)=A ₂ d ² +A ₁ d+A ₀=0 means for obtaining two solutions for the distance to fault according to $\begin{matrix} {d_{1} = \frac{{- A_{1}} - \sqrt{A_{1}^{2} - {4A_{2}A_{0}}}}{2A_{2}}} \\ {d_{2} = \frac{{- A_{1}} + \sqrt{A_{1}^{2} - {4A_{2}A_{0}}}}{2A_{2}}} \end{matrix}$ means for comparing the solutions as 0<(d ₁ or d ₂)<1pu where pu is the length of transmission lines between the line terminals.
 5. Use of a device according to claim 4 to determine the distance to fault in a single transmission line.
 6. Use of a device according to claim 4 to determine the distance to fault in parallel mutually coupled transmission line. 